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Abstract 


The tt charge asymmetry is measured in proton-proton collisions at a centre-of-mass 
energy of 8 TeV. The data, collected with the CMS experiment at the LHC, corre¬ 
spond to an integrated luminosity of 19.7fb^^ Selected events contain an electron 
or a muon and four or more jets, where at least one jet is identified as originat¬ 
ing from b-quark hadronization. The inclusive charge asymmetry is found to be 
0.0010 ± 0.0068 (stat) ± 0.0037 (syst). In addition, differential charge asymmetries as 
a function of rapidity, transverse momentum, and invariant mass of the tf system 
are studied. For the first time at the LHC, the measurements are also performed in 
a reduced fiducial phase space of top quark pair production, with an integrated re¬ 
sult of —0.0035 ± 0.0072 (stat) ± 0.0031 (syst). All measurements are consistent within 
two standard deviations with zero asymmetry as well as with the predictions of the 
standard model. 
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1 Introduction 


The top quark offers an excellent opportunity to search for deviations from the standard model 
(SM), as its large mass makes it unique among all quarks. A possible hint for new physics in 
the top quark sector is the discrepancy between the measured tf forward-backward asymmetry 
and the SM expectations, reported by the CDF |[Tl|2l and DO DsHSl collaborations at the Tevatron. 
Although this discrepancy has become smaller as the measurements and SM calculations ||6HZl 
have been refined, it has generated a number of theoretical explanations invoking contributions 
from physics beyond the SM (BSM). These have in turn led to models based on axigluons or 7! 
bosons as mediators in the tf production process. An overview of the theoretical explanations 
can be found in Ref. I'S] and references therein. 


At hadron colliders top quark pairs are produced predominantly in the processes of gluon- 
gluon fusion and quark-antiquark annihilation. At leading order (LO), the tf production is 
symmetric with respect to the exchange of the top quark and antiquark. At higher orders, QCD 
radiative corrections to the qq —> tf process induce an asymmetry in the differential distribu¬ 
tions of top quarks and antiquarks. The interference between initial- and final-state radiation 
(ISR and FSR) processes, as well as the interference between the Born and box diagrams, gen¬ 
erate a correlation between the direction of the top quark momentum and that of the incoming 
quark |9j. Similarly, the direction of the top antiquark momentum is related to that of the in¬ 
coming antiquark. These processes induce a forward-backward asymmetry (Afb) at the Teva¬ 
tron pp collider. The charge-symmetric pp collisions at the CERN LHC result in a different 
effect. At the LHC, the larger average momentum fraction of the valence quarks leads to an ex¬ 
cess of top quarks produced in the forward and backward directions, while the top antiquarks 
are produced more centrally. This makes the difference in the absolute values of the rapiditie^ 
of the top quark and antiquark, A |j/| = |i/t| — \yi\,^ suitable observable to measure the tf charge 
asymmetry at the LHC experiments. Using the sensitive variable, the charge asymmetry can 
be defined as 


_N+ -N- 
^ ~ N+ + N-' 


( 1 ) 


where N+ and represent the number of events with positive and negative values of A|y|, 
respectively. Theoretical predictions for this observable are of order 1% in the SM IfTOl [TT|, 
but its sensitivity to new physics makes measurements of the effect interesting even when the 
precision is not high enough to establish the existence of the SM charge asymmetry. Both the 
CMS and ATLAS collaborations have published results based on the data collected at a centre- 
of-mass energy = 7TeV, which are in agreement with the SM predictions Ill2l - [T^ . 


To shed light on the possible existence and the nature of new physics contributions, it is crucial 
to measure not only the inclusive asymmetry but also Ac as a function of variables magnifying 
the tf charge asymmetry. For this purpose Eq. Q is modified to consider only events in a 
specific bin of the given variable. 


In this letter, we present an inclusive measurement and three differential measurements of the 
tf charge asymmetry. The three differential variables, which are each sensitive to a different 
contribution to the charge asymmetry, include the tf system rapidity |yjt|, its transverse mo¬ 
mentum pj, and its invariant mass The measurements use the data collected with the CMS 
experiment in 2012 corresponding to an integrated luminosity of 19.7fb^^ at ^/s = 8 TeV. 

The variable \y^^\ is sensitive to the ratio of the contributions from the qq and gg initial states 

^The rapidity is defined as y = (1/2) ln[(£ -|- pz)/(E — pz)], where £ denotes the particle energy and pz its 
momentum component along the counterclockwise beam direction. 
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3 Simulated samples 


to tt production. The charge-symmetric gluon fusion process is dominant in the central region, 
while tf production through qq annihilation mostly produces events with the tf pair at larger 
rapidities, which implies an enhancement of the charge asymmetry with increasing \y^i\ IllOl . 

The ratio of the positive and negative contributions to the overall asymmetry depends on pj. In 
the SM the interference between the Born and the box diagrams leads to a positive contribution, 
while the interference between ISR and FSR results in a negative contribution. The presence 
of additional hard radiation implies, on average, a higher transverse momentum (px) of the tt 
system. Consequently, in events with large values of pj, the negative contribution from the 
ISR-FSR interference is enhanced KTOI . 

The charge asymmetry is expected to depend on since the contribution of the qq initial state 
process is enhanced for larger values of this variable. It is also sensitive to BSM contributions; 
new heavy particles could be exchanged between initial quarks and antiquarks and contribute 
to the tt production (see, e.g. Ref. 1161 and references therein). The amplitudes associated with 
these new contributions would interfere with those of the SM processes, and depending on the 
model they could lead to an increasing tt charge asymmetry with increasing 

Because only a part of the tt phase space is experimentally accessible, measurements of the 
charge asymmetry that are to be compared to theoretical predictions necessarily include an ex¬ 
trapolation to a more well-defined phase space. To this end a fiducial phase space is defined 
that emulates the restrictions of the measurable phase space while allowing for the calculation 
of theoretical predictions. This minimizes the need for extrapolation, which can be subject to 
unpredictable uncertainties if there are significant BSM contributions. An additional extrapo¬ 
lation to the full phase space of top quark pair production is provided as well, which allows 
for an easier comparison to the results of other measurements and theoretical calculations. 


2 CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diam¬ 
eter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and 
strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scin¬ 
tillator hadron calorimeter, each composed of a barrel and two endcap sections. The inner 
tracker measures trajectories of charged particles within the pseudorapidity range \rj\ < 2.5, 
while the calorimeters provide coverage up to |?/| = 3.0. The pseudorapidity is defined as 
rj = — ln(tan0/2), with the polar angle 9 being measured relative to the counterclockwise 
beam direction. The ECAL has an energy resolution of 3% or better for the range of electron 
energies relevant for this analysis. Extensive forward calorimetry complements the coverage 
provided by the barrel and endcap detectors. Muons are measured in the pseudorapidity range 
\i]\ < 2.4 using gas-ionization detectors embedded in the steel flux-return yoke outside the 
solenoid. Matching muons to tracks measured in the silicon tracker results in a relative pj res¬ 
olution for muons with 20 < px < 100 GeV of 1.3-2.0% in the barrel and better than 6% in the 
endcaps. The pj resolution in the barrel is better than 10% for muons with pj up to 1 TeV lllTl . 
A more detailed description of the CMS detector, together with a definition of the coordinate 
system used and the relevant kinematic variables, can be found in Ref. HT^ . 


3 Simulated samples 

For several steps of the measurement, samples of simulated events are used to model both 
the signal process and the background contributions arising from the production of single top 
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quarks and vector bosons in association with jets (W+jets and Z+jets). An additional back¬ 
ground contribution from QCD multijet events is modelled using a template derived from data; 
see Section Top quark pairs are produced with the next-to-leading order (NLO) generator 
POWHEG, version 1.0 IIT9] - t22l . using the CTIO | 23l parton distribution functions (PDF). The 
electroweak production of single top quarks, in the f-channel and in association with a W bo¬ 
son (tW-channel), is simulated using POWHEG and the CTEQ6M PDF set Il24ll . The production 
of electroweak vector bosons in association with jets is simulated using MadGraph, version 
5.1.3.30 |25l, and the CTEQ6L1 El PDF set. 

For the simulation of tf and single top quark events the top quark mass is set to 172.5 GeV. 
For all samples, PYTHIA, version 6.426 1261 , is used for the description of the parton showering 
and hadronization. The simulations include additional proton-proton interactions in the same 
bunch crossing (in-time pileup) and in earlier/later bunch crossings (out-of-time pileup) with 
the same frequency of occurrence as observed in the data. 

Differential cross section measurements Il27ll have shown that the px spectrum of the top quarks 
in tf events is significantly softer than the one generated by the used simulation programs. To 
correct for this effect, the simulated tf sample is reweighted according to scale factors derived 
from these measurements. 


4 Event selection 

The analysis uses tf events in which one of the W bosons from a top quark decay subsequently 
decays into an electron or muon and the corresponding neutrino, and the other W boson de¬ 
cays into a pair of quarks. We therefore select events containing one electron or muon and four 
or more jets, at least one of which is identified as originating from the hadronization of a bot¬ 
tom quark. To be considered for the offline analysis, the events must pass a single-electron or 
a single-muon trigger with pj thresholds of 27 and 24 GeV for the electron and muon, respec¬ 
tively. 

The particle-flow (PF) algorithm ||28j|29| is used to reconstruct electrons, muons, and jets in the 
event. The algorithm reconstructs and identifies each individual particle with an optimized 
combination of information from the various elements of the GMS detector. The reconstructed 
PF candidates are divided into five classes: electrons, muons, photons, charged hadrons, and 
neutral hadrons. 

The primary vertex of the event Il30l is identified as the reconstructed vertex with the highest 
sum of squared transverse momenta of the associated charged particles. For an event to be 
accepted, the primary vertex must satisfy criteria on its location within the detector volume, as 
well as on the quality of its reconstruction. 

Electron candidates are required to have a transverse momentum larger than 30 GeV and be 
within \tj\ < 2.5, excluding the transition region between the EG AT barrel and endcaps of 
1.44 < IpscI < 1-57 since the reconstruction of an electron object in this region is not optimal, 
where rjsc is the pseudorapidity of the electron candidate supercluster ISTl . Furthermore, elec¬ 
tron candidates are selected based on the value of a multivariate discriminant, which combines 
different variables related to calorimetry and tracking parameters, but also px and t] of the 
electron candidate. The electron definition also encompasses a conversion rejection method 
aimed at identifying electrons from photon conversions. Detailed information on the electron 
reconstruction can be found in Ref. 1311 . 

Muons are required to have \ f]\ <2.1 and px > 26 GeV, with further requirements on the qual- 
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5 Definition of a fiduciai phase space 


ity of the muon reconstruction and the compatibility with the primary vertex of the event. The 
rj requirement reflects the coverage of the single-muon trigger. Details on the muon reconstruc¬ 
tion can be found in Ref. Ezl. 

Additionally, electron and muon candidates must be isolated. The isolation is quantified by 
the variable defined as the sum of reconstructed transverse momenta of nearby PF objects 
divided by the lepton transverse momentum (pj), corrected for pileup effects IMI using the 
effective area (in p-cp space) of the lepton and the energy density in the event. Electrons and 
muons are required to have < 0.1 and < 0.12, respectively, using isolation cones with radii 
of 0.3 and 0.4 in t]-(p space. 

Events with additional electrons and muons are vetoed. The lepton veto is based on a looser 
definition of the lepton identification. In this definition, electrons must have px > 20 GeV, 
Ip I < 2.5 and < 0.15, while passing a loose criterion on the value of the multivariate 
discriminant. Muons are required to have px > 10 GeV, |p| < 2.5, and < 0.2. 

Jets are clustered from PE particles with the anti-Zcx 1321 algorithm with a distance parameter 
of 0.5. Gharged hadrons identified as originating from pileup vertices are removed before clus¬ 
tering into jets. Eurther corrections 1331 to the jet energy are applied, depending on jet rj and 
Px, the jet area in t]-(p space, and the median px density of the event. The selected jets must 
lie within |p| <2.5 and are required to have px > 30GeV. The jet px resolution in data is ap¬ 
proximately 10% worse compared to simulations. To account for this, the transverse momenta 
of jets in the simulated samples are smeared accordingly. Einally, jets from the hadronization 
of b quarks are identified using the medium working point of the combined secondary vertex 
algorithm l34l . The b tag identification efficiency of this algorithm for b jets with px > 30 GeV 
and I p I < 2.4 varies between 60 and 70%, while the misidentification rate for jets arising from 
light quarks or gluons is about 1% 1351 . 

With the applied event selection we find a total of 171121 events with an electron in the final 
state, hereafter referred to as the electron+jets channel, and 192123 events in the muon+jets 
channel. 


5 Definition of a fiducial phase space 

Because of the offline event selection, only a subset of the events collected by the GMS detec¬ 
tor is used in the analysis. To allow for a comparison of the measurements with the theoretical 
calculations, an extrapolation to a well-defined phase space needs to be performed. The extrap¬ 
olation relies on a correct modeling of the ratio of the number of events in the measured phase 
space to that in the extrapolated one; such a ratio, however, may be affected by new physics. 
The simple approach, which is extrapolation to the full phase space of tf production, entails a 
large dependence on the model assumptions. 

As an alternative, a fiducial phase space is defined using generator-level selection criteria that 
mimic the reconstruction-level criteria applied during the nominal selection. The ratio of the 
number of fiducial events to the number of reconstruction-level selected events, determined 
from simulation, is then applied to the data to estimate the distribution of an observable in the 
fiducial region. 

Because of the physical and topological similarity of events in the selected and fiducial phase 
spaces, new physics contributions are expected to affect both in approximately the same way, 
leaving the ratio unchanged. Thus this extrapolation to the fiducial phase space is nearly 
model-independent. It should be noted that the desired model-independence is achieved by 
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using a purely multiplicative correction; a prior subtraction of non-fiducial tf events in the se¬ 
lected phase space would require a larger reliance on the model assumptions, as there would 
be no cancellation of uncertainties. 

Jets of generated particles in simulated events are used to emulate the selection steps acting 
on reconstructed jets. Hadron-level particles are clustered into jets using the anti-Zcj algorithm 
with a distance parameter of 0.5, as used for the reconstructed jets. The clustering includes 
charged leptons and neutrinos, except those created in the leptonic decay of W bosons origi¬ 
nating from top quarks. It should be noted that the selection criteria for charged leptons are 
applied only to leptons originating from top quark decays. 

Using these objects the following selection requirements are applied. The event needs to con¬ 
tain exactly one electron (or muon) with pj > 30 (26) GeV and \tj\ < 2.5 (2.1). Any event that 
contains an additional electron (or muon) with pj > 20 (10) GeV and \fi\ < 2.5 is rejected. 
At least four generator-level jets with pj > 30GeV, \rj\ < 2.5 are required. The event is re¬ 
jected if the axes of any such jets have an angular separation of AR < 0.4 to the lepton, where 
AR = V(A/;)2 + (A(J))2 is calculated using the differences in the azimuthal angles (p and pseu¬ 
dorapidities p. This criterion serves as an emulation of the lepton isolation criteria, which use 
similar radii and are hard to implement for theoretical calculations. 

The fiducial region contains about 10% of the events of the full phase space. Roughly 50% 
of the events in the fiducial region pass the selection outlined in Section with the largest 
inefficiencies caused by the lepton selection and trigger requirements. In comparison, only 
1.5% of the events outside the fiducial region fulfil the event selection criteria, making up about 
20 % of the selected events. 


6 Estimation of background contributions 

For the estimation of the background contributions we make use of the discriminating power 
of the transverse mass of the W boson, , and of M 3 , the invariant mass of the combination of 
three jets that corresponds to the largest vectorially summed pj. This estimation is necessary for 
the subtraction of the background contributions of the measured data, as described in Section]^ 
The variable is calculated from the transverse momentum of the charged lepton pj and the 
missing transverse momentum vector The latter is defined as the pj imbalance of the 

reconstructed PF objects, taking into account the propagation of jet energy corrections to this 
observable. Its magnitude is referred to as Neglecting the lepton masses, is defined 

as 

• p^j). ( 2 ) 

The background estimation is made with a binned maximum likelihood fit for each lepton 
channel. A simultaneous fit in and M 3 is performed in two disjoint data samples, corre¬ 
sponding to < 50 GeV and >50 GeV. The distribution is fitted in the low-m™ sample, 
which is rich in QGD multijet events and yields a good discrimination between the QGD mul¬ 
tijet process and processes containing a genuine W boson. The distribution of M 3 is not as 
dependent on the choice of event sample; it is fitted in the complementary high-iw^ sample to 
avoid using the same events for both fits. 

For the tf, W+jets, Z+jets, and single top quark processes, simulated samples are used to model 
the shapes of the ni^ and M 3 distributions. The contribution from multijet background events 
is estimated from data control samples containing nonisolated or poorly identified leptons. 
Rate constraints corresponding to Gaussian uncertainties of 20% are introduced into the like- 
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7 Measurement of the tt charge asymmetry 


along with 
been taken 
The uncer- 
constraints 


Process 

Electron+jets 

Muon+jets 

Single top quark (f + tW) 

7016 ± 1328 

7302 ± 1663 

W+jets 

22 508 ±1460 

20 522 ±1606 

Z+jets 

2345 ± 510 

2046 ± 415 

QCD multijet 

6136 ± 1201 

4199 ± 588 

Total background 

38005 ±1491 

34096 ±1495 

tt 

133130 ± 1521 

158 058 ±1538 

Observed data 

171121 

192123 


Table 1: Number of events for background and tt contributions from fits to data, 
their statistical uncertainties. The correlations between the individual values have 
into account for the determination of the uncertainty on the total background yield, 
tarnties quoted for the single top quark and Z+jets backgrounds are driven by the 
used as inputs for the likelihood fit. 
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Figure 1: Comparison of the combined lepton+jets data with simulated contributions for the 
distributions in ni^ and M 3 . The simulated signal and background contributions are normal¬ 
ized to the results of the fits in Table The last bin in each histogram includes the overflow 
values. Additionally, the ratio of the data to the sum of the simulated contributions is shown, 
with the statistical uncertainties of the simulated contributions (including the uncertainties in 
the fit) indicated by the blue hatched region. 


lihood function for the Z+jets and single top quark processes according to the respective NLO 
cross sections, while the rates of the other processes are free parameters of the fit. The width of 
the constraints is motivated by the uncertainties of measurements and theoretical calculations 
of the corresponding cross sections Il36ti40l . A detailed description of the fitting procedure can 
be found in Ref. |4T] . 

Table [^summarizes the results of the fits. Figurej^shows the two variables used for the estima¬ 
tion of the background, with the individual simulated contributions normalized to the results 
from the fit. 


7 Measurement of the tt charge asymmetry 

The measurement of the tt charge asymmetry is based on the reconstructed four-momenta of 
the t and f quarks in each event. We reconstruct the leptonically decaying W boson from pj 
and pjF'®® and associate the measured jets in the event with quarks in the tt decay chain. The 
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association is done using a likelihood criterion based on the b tagging discriminator values 
of the jets and the corresponding reconstructed masses of the top quarks and W bosons. The 
reconstruction procedure is described in detail in Ref. |4T|. 


The reconstructed top quark and antiquark four-momenta are used to obtain the inclusive and 
differential distributions of A|y|, and the charge asymmetry is calculated from the number of 
entries with A|i/| > 0 and < 0. In case of the differential measurements, the asymmetries are 
calculated separately for the different bins in the kinematic variable Vi, where Vi is either |i/fj|, 
or OTtf. 

To allow for a comparison of the resulting asymmetry and the predictions from theory, the 
reconstructed distributions of A|i/| and the three kinematic variables are corrected for back¬ 
ground contributions, resolution, and selection efficiency. 

In the first correction step, the distributions of background processes, as used in Section]^ are 
normalized to the estimated rates (see Table and subtracted from the data, assuming Gaus¬ 
sian uncertainties in the background rates as well as in statistical fluctuations in the background 
templates. The correlations among the individual background rates are taken into account. 
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Figure 2: Migration matrix between generated (A|y|gen) and reconstructed (A|y|rec) rapidity 
differences (top left) and selection efficiency with respect to the full phase space as a function 
of AI j/1 gen (top right) of the inclusive measurement. Migration matrix (bottom left) and selection 
efficiency (bottom right) for the measurement differential in m^l. 


The resulting background-subtracted distributions are translated from the reconstruction level 
to parton level within the phase space of the selected events. Afterwards, acceptance correc¬ 
tions are applied, correcting either to the fiducial phase space described in Section or to the 
full phase space. Apart from this last step, the measurements for both phase spaces are identi- 
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7 Measurement of the tt charge asymmetry 


cal. After the corrections have been applied, the resulting distributions are independent of the 
detector and analysis specifications. 

The above corrections are obtained by applying an unfolding procedure to the data Il42l through 
a generalized matrix inversion method. In this method, the resolution and selection effects are 
described by a response matrix R that translates the true spectrum x into the measured spec¬ 
trum w = Rx. As reconstruction and selection effects factorize, the response matrix R can be 
seen as the product of a migration matrix, describing reconstruction effects, and a diagonal 
matrix containing the selection efficiencies, describing acceptance effects. Both the migration 
matrix and the selection efficiencies are determined from simulated tt events. As the compo¬ 
nents corresponding to the electron+jets and muon+jets channels are found to be very similar, 
they are combined to yield a method that can be applied to the summed data of both chan¬ 
nels. In this combination the individual components are scaled according to the scale factors 
obtained via the background estimation. The unfolding procedure used in the inclusive mea¬ 
surement, described in detail in Ref. 1411 . is also used for the two-dimensional distributions of 
the differential measurements. 

This analysis uses 12 bins for the unfolded A|i/| distribution in the inclusive measurement and 
8 bins for the same distribution in the differential measurements. The unfolded Vj distributions 
use 3 bins, with one additional measurement being performed using 6 bins in The addi¬ 
tional measurement provides finely grained results in the region of high The ranges for the 
bins in these distributions are given in Table It should be noted that the outermost bins of 
A|y| extend to infinity. 

Table 2: The bin ranges for the individual bins of the differential measurements. Two different 
choices of binning are used for the distribution of m^i. 


Bin 

ll/tti 

pf (GeV) 

mjt (GeV) 

nin (GeV) 

1 

0-0.34 

0-41 

0-430 

0-420 

2 

0.34-0.75 

41-92 

430-530 

420-500 

3 

0.75-00 

92-00 

530-OO 

500-600 

4 




600-750 

5 




750-900 

6 




900-OO 


In the corresponding reconstructed spectra the numbers of bins along both axes are doubled, 
as is recommended for the applied unfolding procedure 1421 . The choice of the bin edges for 
A|y| is different in each bin of Vi, resulting in different amounts of vertical overlap between 
horizontally neighbouring bins in the two-dimensional distributions (for illustration see the 
binning in Fig. bottom right). 

To limit the magnification of statistical uncertainties due to the unfolding procedure, a regular¬ 
ization is applied that suppresses solutions with large fluctuations between neighbouring bins. 
The strength of the regularization is determined by minimizing the statistical correlations be¬ 
tween bins in the unfolded spectrum. Different strengths are used for the regularization along 
the sensitive variable within each bin of the kinematic variable. Similarly, the regularization 
along the kinematic variable is adjusted separately for each bin of the kinematic variable. 

Separate migration matrices are used for the inclusive measurement and for each of the dif¬ 
ferential measurements. Figure]^ shows the migration matrices for the inclusive measurement 
and, as an example, for the differential measurement in m^i. For the inclusive measurement 
the migration matrix describes the migration of selected events from true values of A|i/| to the 
reconstructed values. For the migration matrices of the differential measurements not only the 
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migration between bins of A|i/| has to be taken into account, but also the migration between 
bins of Vj. For a measurement in 3 unfolded bins of Vi these migration matrices feature a grid 
of 6 X 3 bins in Vi, with each of these bins representing a 16 x 8 migration matrix describing the 
migration between different A|i/| values. 

The values of A|i/| and h, also affect the probability for an event to fulfil the event selection 
criteria. The selection efficiencies relative to the full phase space for the inclusive measurement 
and for the differential measurement in are depicted in Fig. The selection efficiency of 
the fiducial phase space is defined by the ratio of all selected events to the events present in 
the fiducial phase space. It should be noted that the selected events also include events that do 
not pass the criteria of the fiducial phase space; their influence is implicitly corrected for in the 
acceptance correction because of the way the selection efficiency is defined. Thus this correction 
is multiplicative in nature, which is justified by the inherent similarity of these events to the 
events that are intended to be measured. 

One limiting factor for the precision of the analysis is the presence of sizeable statistical fluc¬ 
tuations in the response matrices as they are obtained from simulated events. To mitigate this 
effect, one can exploit an approximate symmetry of the response matrix under charge conjuga¬ 
tion. For events resulting from a charge-symmetric initial state like gluon-gluon fusion it can 
be assumed that reconstruction effects also have a predominantly charge-symmetric behaviour. 
From this reasoning, the symmetry is enforced for this analysis by averaging those bins of the 
gluon-gluon contribution to the response matrix that correspond to each other under charge 
conjugation. 

The correctness of the unfolding procedure has been verified with pseudo-experiments, each 
of which provides a randomly generated sample distribution from the templates used in the 
analysis. 


8 Estimation of systematic uncertainties 

The measured charge asymmetry Aq is affected by several sources of systematic uncertainty 
Effects altering the direction of the reconstructed top quark momenta can change the value 
of the reconstructed charge asymmetry. Systematic uncertainties with an impact on the dif¬ 
ferential selection efficiency, as well as variations in the rates and modelling of background 
contributions, can also bias the result. To evaluate each source of systematic uncertainty, a new 
background estimation is performed and the measurement is repeated on data using modified 
simulated samples. The differences in unfolded asymmetries are then used to construct a sys¬ 
tematic asymmetry covariance matrix in a loose analogy to statistical covariance matrices. For 
an uncertainty described by a single systematic shift a covariance of 

COv(x, y) = {x- Xnom) (}/ “ 1/nom) (3) 

is used, with x and y referring to bins of the asymmetry distribution resulting from the system¬ 
atic shift and Vnom and r/nom being the results of the nominal measurement. For uncertainties 
that are determined using exactly two variations (indexed by 1 and 2) the absolute values of 
the maximal shifts observed in each result bin, Ax^ax and Ai/max/ are determined separately; 
the covariance is then defined as 

cov(x,y) = AxmaxAymax sign((xi -X2) (yi -l/l))- (4) 

This procedure corresponds to a symmetrization of the largest observed shifts and thus con¬ 
stitutes a more conservative uncertainty estimate than an approach based on a direct analogy 
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8 Estimation of systematic uncertainties 


with statistical covariance definitions. The covariance matrices of all systematic uncertainties 
are added up to yield a resultant matrix where the diagonal elements are the variances. 

In the following, a summary of the studied sources of systematic uncertainty is given. 

The corrections to the jet energy scale and jet energy resolution are varied within their tj- and 
pj-dependent uncertainties to estimate their effects on the measurement. The effect of varia¬ 
tions in the frequency of occurrence of pileup events is determined using reweighted simulated 
samples. Differences between data and simulations in the b tagging efficiency and the lepton 
selection efficiency are determined as scale factors that depend on pj, or p and respectively. 
Effects due to uncertainties on these scale factors are studied by varying them as a function of p 
within their uncertainties and, in the case of the lepton selection efficiency, also as a function of 
the lepton charge. The effect of lepton charge misidentification is very small and is neglected. 

To estimate the influence of a possible mismodelling of the simulated W+jets background, 
the measurement is repeated using a W+jets template determined from a sideband region in 
data, defined by an inversion of the requirement of a b-tagged selected jet. The template is 
reweighted to account for the differences between the signal and sideband regions, which are 
determined from the simulation. 

The uncertainty in the multijet background modelling in the electron+jets channel is deter¬ 
mined by replacing the nominal template, which is estimated using two sideband regions de¬ 
fined either by inverted isolation or by inverted identification criteria, with templates derived 
from only one of the sideband regions each. Meanwhile, in the muon+jets channel, only the 
template from the isolation-inverted sample can be used, so a conservative estimation of the 
uncertainty in this background contribution is performed by taking the maximum deviation 
out of three scenarios where the multijet template is replaced with the tt signal template, with 
the simulated W+jets template, or with a template obtained by inverting the sign of the sensi¬ 
tive variable in the multijet template itself. 

In contrast to the other systematic effects, the uncertainty due to the unfolding method is esti¬ 
mated by unfolding simulated samples instead of data. The simulated tt events are reweighted 
to reproduce the observed asymmetries in the differential measurements based on data, and 
the resulting reconstruction-level spectra are unfolded. The deviations between the unfolded 
asymmetries and the reweighted true asymmetries are taken to be a measure of the model 
dependence of the unfolding procedure in the observed point in phase space. The actual un¬ 
certainty of each measurement is estimated as the square root of the average squared deviations 
produced by the unfolding in the three reweighting scenarios corresponding to the three kine¬ 
matic variables. 

To estimate the uncertainty resulting from possible mismodelling of the tt signal, samples 
of simulated tt events produced with MadGraph are compared to samples produced with 
POWHEG, both interfaced to PYTHIA for the modelling of the parton shower. In a similar way 
the impact of a possible mismodelling of parton shower and hadronization is studied by us¬ 
ing HERWIG Il43ll44l . as opposed to PYTHIA, for the simulation of the signal, with the hard- 
scattering matrix element being simulated by either POWHEG or MC@NLO Il45l . As a measure 
of the uncertainty related to the performed reweighting as a function of the top quark pj, de¬ 
scribed in Section]^ the measurement is repeated using samples without reweighting. Finally, 
the impact of variations in the renormalization and factorization scales (/Ir and /ip) in the simu¬ 
lated tt events is determined using dedicated samples generated at scales varied up and down 
by factors of 2. 

The systematic uncertainty on the measured asymmetry from the choice of PDFs for the collid- 
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ing protons is estimated using the LHAPDF package and the uncertainty in the CTIO Il23l . 
MSTW2008 BZl, and NNPDF2.1 ESl PDF sets. 

The contributions of the different sources of systematic uncertainties to the total uncertainty 
of the inclusive measurements are summarized in Table The table also shows the ranges of 
systematic uncertainties in the 3-binned differential measurements to illustrate the magnitudes 
of the individual contributions. Because the measurements in the two phase spaces differ only 
by the acceptance corrections, the uncertainties can be seen to behave similarly for the two 
cases. 

Table 3: Uncertainties for the inclusive measurement of Aq and ranges of uncertainties for the 
differential measurements using three bins for the kinematic variable. Numbers are given for 
measurements in the fiducial phase space (fid. PS) and in the full phase space (full PS). 


Uncertainty source 

Inclusive Aq uncertainty Differential Aq uncertainty 
fid. PS full PS fid. PS full PS 

Jet energy scale 

0.0020 

0.0018 

0.0009-0.0066 

0.0008-0.0063 

Jet energy resolution 

0.0003 

0.0003 

0.0005-0.0020 

0.0005-0.0020 

Pileup 

0.0006 

0.0006 

0.0002-0.0027 

0.0003-0.0027 

b tagging 

0.0009 

0.0008 

0.0002-0.0033 

0.0002-0.0032 

Lepton selection efficiency 

0.0009 

0.0009 

0.0005-0.0016 

0.0005-0.0017 

W+jets background 

0.0005 

0.0007 

0.0003-0.0030 

0.0005-0.0025 

QCD multijet background 

0.0010 

0.0009 

0.0008-0.0030 

0.0011-0.0028 

Unfolding 

0.0012 

0.0022 

0.0004-0.0023 

0.0011-0.0033 

Generator 

0.0002 

0.0005 

0.0008-0.0058 

0.0007-0.0043 

Hadronization 

0.0010 

0.0011 

0.0007-0.0046 

0.0008-0.0040 

Top quark pj reweightmg 

0.0000 

0.0002 

0.0000-0.0014 

0.0001-0.0015 

Pr and f/p scales 

0.0002 

0.0007 

0.0008-0.0057 

0.0009-0.0064 

PDF 

0.0002 

0.0003 

0.0004-0.0014 

0.0004-0.0012 

Total syst. uncertainty 

0.0031 

0.0037 

0.0043-0.0120 

0.0041-0.0115 

Statistical uncertainty 

0.0072 

0.0068 

0.0078-0.0181 

0.0078-0.0172 

Total uncertainty 

0.0078 

0.0077 

0.0094-0.0217 

0.0094-0.0207 


9 Results 

Table l^gives the values of the measured inclusive asymmetry at the different stages of the anal¬ 
ysis, while the unfolded A|i/| distributions for the fiducial and full phase spaces are shown in 
Kg-0 The latter two distributions are shown in the form of normalized differential cross sec¬ 
tions as a function of A |y |. All unfolded quantities correspond to the parton level. The statistical 
uncertainty of all quoted results encompasses the subdominant effects of the limited number of 
simulated events used for the measurement. It should be noted that the acceptance corrections 
for the two phase spaces differ as a function of A|i/|; as a result, statistical fluctuations in the 
data can have different effects on the measured asymmetries. 

The uncertainty in the theoretical prediction by Kiihn and Rodrigo lllOl is estimated by varying 
the top quark mass, the PDFs, and the f/R and f/p scales, with the scale uncertainties being the 
dominant effect. The uncertainty in the theoretical prediction by Bernreuther and Si IfTTl |4^ 
contains only the effects of variations of the and scales. The tf charge asymmetry for the 
fiducial phase space is computed with the tt production and semileptonic/non-leptonic tt de¬ 
cay matrix elements at NLO. The top quark decay matrix elements at NLO contain additional 
scale dependencies. This results in a larger scale uncertainty as compared to the charge asym- 
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9 Results 


metry for the full phase space. Another recent CMS analysis of the inclusive charge asymmetry 
in the full phase space IISOll . which uses a slightly more model-dependent approach to achieve 
lower uncertainties, and a recently published ATLAS measurement of inclusive and differential 
charge asymmetries USTII both yield results that are comparable to the ones presented here. 


Table 4: The measured inclusive asymmetry at the different stages of the analysis and the 
corresponding theoretical predictions from the SM. 


Asymmetry (Aq) 

Reconstructed 

Background-subtracted 

Corrected for migration effects 

0.0036 ± 0.0017 (stat) 
0.0008 ± 0.0023 (stat) 
-0.0042 ± 0.0072 (stat) 

Fiducial phase space 

Theoretical prediction [Bernreuther, Si] liTTllT^ 

-0.0035 ± 0.0072 (stat) ± 0.0031 (syst) 
0.0101 ±0.0010 

Full phase space 

Theoretical prediction [Kiihn, Rodrigo] KTOl 
Theoretical prediction [Bernreuther, Si] IITTll4^ 

0.0010 ± 0.0068 (stat) ± 0.0037 (syst) 
0.0102 ± 0.0005 

0.0111 ±0.0004 

19.7 fb ' (8 TeV) 

19.7fb ' (8 TeV) 


g 0.8: 

T3 
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Figure 3: Unfolded inclusive A|i/| distribution in the fiducial phase space (left) and in the full 
phase space (right). The uncertainties on the data points represent the statistical uncertainties 
due to the limited amounts of data and simulated events. The measured values are compared 
to NLO predictions for the SM based on calculations by Kiihn and Rodrigo (K&R) llTOl and 
Bernreuther and Si (B&S) 1(11114^ . 


The results of the differential measurements in the fiducial phase space are shown in Fig. 
and the extrapolation to the full phase space in Fig. The measured values are compared to 
predictions from SM calculations IITOl ITTl |4^ as well as to predictions from an effective field 
theory Il52l |53]|. The latter is capable of reproducing the CDF results |2| by introducing an 
anomalous effective axial-vector coupling to the gluon at the one-loop level. The gluon-quark 
vertex is treated in the approximation of an effective field theory with a scale for new physics 
contributions of order 1.5-2.0 TeV. Predictions for the asymmetry as a function of pj are not 
available for this theory and for one of the SM calculations. Because of the importance of the 
region of high m^^ for the detection of new physics, we provide an additional, more finely- 
grained differential measurement of the asymmetry as a function of this observable. 

Both of the inclusive results yield values that are slightly smaller than the SM predictions, with 
the larger deviation being in the fiducial phase space and corresponding to 1.7 standard devia¬ 
tions. The differential measurements show a good agreement with the SM predictions. For the 





















13 


o 

< 


19.7 fb ' (8 TeV) 


0.5 


|y. 


50 


o 

< 


0 . 02 - 

0 - 

- 0 . 02 - 

-0.04- 


19.7fb ' (8 TeV) 



100 150 

p« (GeV) 


_CWIS 

- l-i-jets, fiducial PS 

1 1 1 1 1 

t Data 

QCD (NLO), B&S 

o 0.1 
< 

CMS 

l+jets, fiducial PS 

- 

II 

1 Data 
- QCD (NLO) 

B&S ; 

- 





0.05 

- 






- 

- 






- 






- 

- 





0 

- 






- 

- 

J _,_ 


J_ , _ ^ 

_^ ^ _1_ 

-0.05 

- 


1 



_,_1_ 

- 


400 


600 


800 

m„- (GeV) 


500 


1000 

m, (GeV) 


Figure 4: Corrected asymmetry as a function of \y^i\ (upper left), pj (upper right), and 
(lower left and lower right). The latter is shown in two different binnings. All results corre¬ 
spond to the fiducial phase space. The measured values are compared to an NLO prediction for 
the SM based on calculations by Bernreuther and Si (B&S) lUTllT^ . The inner bars indicate the 
statistical uncertainties, while the outer bars represent the statistical and systematic uncertain¬ 
ties added in quadrature. 


benchmark model involving an effective axial-vector coupling of the gluon, the measurement 
at high excludes new physics scales below 1.5 TeV at the 95% confidence level. 


10 Summary 

Inclusive and differential measurements of the charge asymmetry in tf production at the LHC 
are presented. The data sample, collected in proton-proton collisions at = 8 TeV with the 
CMS detector, corresponds to an integrated luminosity of 19.7 fb^^. Events with top quark 
pairs decaying into the electron-i-jets and muon-i-jets channels are selected and a full tf event 
reconstruction is performed to determine the four-momenta of the top quarks and antiquarks. 
The observed distributions are then corrected for acceptance and reconstruction effects. For the 
first time at the LHC, acceptance corrections to the tf charge asymmetry are performed not only 
to the full phase space but also to a fiducial phase space. Within two standard deviations, all 
measured values are consistent with the predictions of the standard model and no hint of new 
physics contributions is observed. The charge asymmetry in the high-mass region is about two 
standard deviations below the predictions from an effective field theory with the scale for new 
physics at 1.5 TeV. 
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10 Summary 





Figure 5: Corrected asymmetry as a function of \y^i\ (upper left), pj (upper right), and 
(lower left and lower right). The latter is shown in two different binnings. All results corre¬ 
spond to the full phase space. The measured values are compared to NLO predictions for the SM 
based on calculations by Kiihn and Rodrigo (K&R) |[T0l and Bernreuther and Si (B&S) lUTl|4^ , 
as well as to the predictions of a model featuring an effective axial-vector coupling of the gluon 
(BAG) Il5^l53l . The inner bars indicate the statistical uncertainties, while the outer bars repre¬ 
sent the statistical and systematic uncertainties added in quadrature. 
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